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ON , Abstract. Most gravitational lens galaxies are early- type galaxies in 

relatively low density environments. We show that they lie on the same 
fundamental plane as early-type galaxies in both local and distant rich 
clusters. Their surface brightness evolution requires a typical star forma- 
Q ■ tion epoch of zj ~ 2-3, almost indistinguishable from that of rich cluster 

Q ' galaxies at comparable redshifts. The restricted galaxy type range of the 

lenses means that photometric redshifts work well even with only 1-3 
filter photometry. We make preliminary measurements of the mass and 
luminosity functions of the lens galaxies, and find they are consistent with 
!>■ ' the standard model used for deriving cosmological limits using lens statis- 

tics. As expected for a mass- weighted sample, they are more massive and 
more luminous than the overall early-type galaxy population. 
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Introduction 



M-f Gravitational lens galaxies are a unique sample, because they are the only galax- 

q ■ ies selected based on mass rather than luminosity. As such, the average prop- 

pi | erties of the lens galaxies are identical to the mass-weighted average properties 

of all galaxies at a given redshift. At intermediate redshifts (0 < Z\ < 1), they 
are also the largest sample of galaxies with known masses outside the cores 
of rich clusters. Historically, the interpretation of the mass measurements has 
k> | been hampered by a lack of accurate measurements of the light. Eliminating 

5_j ■ this rather peculiar problem for an astronomical sample is a primary goal of the 

CASTLES (CfA/Arizona Space Telescope Lens Survey) project (see Falco et al. 
in these proceedings). 

If we combine accurate measurements of the mass and the light with the 
wide range of lens redshifts, then we have an excellent tool for studying the 
evolution of galaxy mass-to-light (M/L) ratios with redshift. Our method is 
similar to that introduced by van Dokkum & Franx (1996) and used by Kelson 
et al. (1997, 1999), van Dokkum et al. (1998, 1999), Pahre et al. (1999ab), 
and Jorgensen et al. (1999) to measure the evolution of early-type galaxies in 
rich clusters over the same redshift range. The lens galaxies are in far lower 
density environments than those in rich clusters - a "field" rather than a cluster 
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population - and theoretical models (e.g. Kauffmann 1996, Kauffmann & Char- 
lot 1998) predict that they should have significantly younger stellar populations 
than galaxies in rich clusters. Studies of the colors and luminosity-effective ra- 
dius correlations (e.g. Schade et al. 1996, 1999, Ziegler et al. 1999), and one 
small fundamental plane sample (Treu et al. 1999) all suggest that there is little 
difference between the field and cluster populations at z ~ 0.5. With the lens 
galaxies we can directly measure the surface brightness evolution to z — 1. 

In this review we summarize our recent results on galaxy evolution (Kochanek 
et al. 1999) and illustrate new results on the mass and luminosity functions of 
the lens galaxies. In §2 we introduce the fundamental plane (FP) of lens galaxies, 
which we use to determine the star formation epoch in §3 and as a photometric 
redshift estimator in §4. In §5 we determine the mass and luminosity function 
of the lenses as compared to local galaxies and simple statistical models, and in 
§6 we review our goals for the future. 

2. The Fundamental Plane And Galaxy Evolution 

The fundamental plane (FP) is a tight correlation between the central velocity 
dispersion a c , intermediate axis effective radius r e , and central surface brightness 
\x e discovered by Djorgovski h, Davis (1987) and Dressier et al. (1987). The FP 
is a consequence of virial equilibrium, the nearly self-similar photometric and 
kinematic properties of early-type galaxies, and systematic rather than random 
variations in the metallicity and stellar population ages with mass or luminosity. 
As a tool for studying galaxy evolution, its key feature is that given any two of 
the three variables one can predict the remaining one. Traditionally this is done 
by comparing the directly measured effective radius r e to that predicted by the 
FP, 

log(rf p //i^ 1 kpc) = 1.24 log (cj c /km s" 1 ) + 0.33^ e (0) - 8.66 (1) 

using Jorgensen, Franx & Kjaergaard (1993, 1995ab, 1996, collectively referred 
to as JFK) renormalized to the closest HST filter (F606W). The only (rapidly?) 
evolving quantity appearing in the FP is the surface brightness, fi e ( z ), through 
the time variation in the average luminosity of a stellar population. Hence, the 
standard FP is normalized using the surface brightness today, [i e {z = 0). 

The key to studying galaxy evolution is the strong variation with redshift 
of iJLe{z) as compared to the invariance of r e and a c . If we observe a galaxy in a 
series of filters j, we can measure the surface brightness /ij, which evolves as 

fij(z) = fij(0) + 10 log (1 + z) + ej (z) + kj(z) (2) 

where /ij(0) is the surface brightness at z = 0, the first redshift term is the 
(1 + z) 4 cosmological dimming, ej(z) is the evolution correction for filter j, and 
kj(z) is the K-correction from the rest frame to the observed frame. If we can 
estimate the "final" surface brightness, /Xj(0), then we can directly measure the 
evolutionary terms. By simply reordering the standard form of the FP, we have 

/if p (0) = -3.76 log(CT c /km s" 1 ) + 3.03 log(r e /h^kpc) + 26.26, (3) 

so if we measure the effective radius and velocity dispersion of an early-type 
galaxy at any redshift, we can predict the surface brightness it will have today 
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Figure 1. The FP of lens galaxies transformed to zero redshift. The 
cosmologies (from left to right) are the Oo = 1.0 flat, Qq = 0.3 open and 
Oo = 0.3 flat models. An instantaneous burst star formation history 
is used with star formation redshifts (from top to bottom) of zj = 1, 
2, 3 and 5. The filled squares are for the lenses with known redshifts 
and for the open squares we have used photometric redshifts (see §4). 
The solid line marks the FP of the local comparison sample; 90% of 
the galaxies in the local JFK sample lie between the dashed lines. The 
mean residual (A = (log(r e /Vg ))) and its dispersion are shown in the 
upper left corner of each panel. These are calculated using only the 
systems with known lens redshifts. 



at z = 0. Thus, for any early- type galaxy we have that 

e, (z) + kj (z) = H (z) - [fif P (0) + 10 log(l + zj\ . (4) 

We can determine ej(z) alone by measuring colors and interpolating to a fixed 
rest wavelength (van Dokkum & Pranx 1996). The uncertainty for any one 
galaxy is not just the measurement uncertainty, as we must include the signifi- 
cant surface brightness scatter of the local FP (0.23 mag/arcsec 2 in JFK). 

Our procedures differ somewhat from the cluster studies because we must 
consider individual lens galaxies rather than samples of cluster galaxies at a 
common redshift. The cluster studies measured the properties of enough galaxies 
in each cluster (6 on average) to use the standard FP (eqn. 1), measure the 
change in the zero-point (the 8.66) with redshift and then interpret it as evolution 
in the population averaged surface brightness or mass-to-light ratio. Since we 
cannot construct directly the FP of lens galaxies at any redshift we are driven 
to our alternate, galaxy by galaxy formulation. We believe, however, that our 
formulation of the analysis is the more physical. The existence of a thin FP 
at higher redshifts is observational evidence that early-type galaxies of similar 
mass have small age differences, but it is not required for the analysis. 

For most lenses we do not have velocity dispersion estimates. We instead 
infer the velocity dispersion from the image separation A6 assuming a singular 
isothermal sphere (SIS) mass distribution. Dynamical models of galaxies in SIS 
halos (Kochanek 1994) demonstrated that the dark matter dispersion, which 
we measure with A8, is coincidentally almost identical to the central stellar 
dispersion. Adopting a normalization scale of a Cif = 225 km s _1 for an L* 
early-type galaxy (as required by models of the distribution of lensed image 
separations, Kochanek 1996, Falco et al. 1998), our velocity dispersion estimate 
becomes 

a c = (225/f)(Ae/2'.'91) 1 / 2 (D s/D LS ) l / 2 km s"\ (5) 

where / = (Jdark/&* — 1-0 ± 0.1 is the normalizing factor between the velocity 
dispersion of the dark matter and the stars. We can also use the FP to estimate 
that its value is / = 1.06 ± 0.07, providing further confirmation of the earlier 
dynamical and lensing results. We explored constant M/L dynamical models as 
well, but they generally fit the data poorly. 

Our final current consists of 29 lenses, 12 of which are missing lens redshifts 
and 6 of which are missing source redshifts. Missing redshifts are the bane of lens- 
ing studies, and we discuss photometric redshift estimates for the lens galaxies 
in §4. The estimates of galaxy evolution and the photometric redshift estimates 
for the lens galaxies are insensitive to our assumptions about unmeasured source 
redshifts. We also analyzed a comparison sample of 54 galaxies from the clus- 
ter FP studies, so that we could compare the evolution of the two populations 
directly. For theoretical models we used the GISSEL96 versions of the Bruzual 
& Chariot (1993) spectral evolution models assuming Hq = 65 km s _1 Mpc -1 . 
For simplicity we have compared the data to solar metallicity models with a 
single burst of star formation at redshift Zf. Models with extended bursts are 
also consistent with the data but require earlier star formation epochs. 
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Figure 2. Evolution and K-corrections for the lens (left) and cluster 
(right) samples in the V=F555W (bottom), I=F814W (middle) and 
H=F160W (top) bands as a function of redshift for a flat Qq = 0.3 
cosmological model. The zero-redshift colors were left in to separate the 
curves. The points are averages in redshift bins with edges at z = 0.25, 
0.50 and 0.75. The error bar on the E+K correction is the standard 
deviation of the galaxies in the bin, not the uncertainty in the mean. 
For the lenses, the filled points are the averages using only the lenses 
with known redshifts, while the open points include all lenses. The 
dashed curves are the no evolution models for each filter, and the solid 
curves are the instantaneous burst models with star formation redshifts 
(from bottom to top) of Zf = 10, 3, 2, 1.5 and 1.0 respectively. 



3. The Star Formation Epoch of Early- Type Galaxies 

The lens galaxies at any given redshift are a mass-weighted sample of all galax- 
ies at that redshift. This means that the mean star formation epoch of the lens 
galaxies is closely related to the mean, mass-weighted star formation epoch of 
all galaxies. Hence, the star formation epoch of the lens galaxies is a far more 
fundamental measurement than the star formation epoch of galaxies in rich clus- 
ters. While we consider only early-type lenses here, they are the vast majority 
of all lenses. Moreover, even the late-type lenses tend to be fairly red and show 
few signs of active star formation. The absence of the numerous, luminous, blue, 
low-mass star forming galaxies in the lens sample confirms that these galaxies 
make a negligible contribution to the total mass in stars. 

We can show that the lens galaxies lie on a coherent FP by evolving their 
properties forward in time and placing them on the present day FP. The trans- 
formations depend on the star formation epoch, Zf, and the cosmological model 
(f^o and Ao). Figure 1 shows the results. The mean logarithmic residual, 
A = (logfe/fg ), is computed only for the lenses with known redshifts. As 
van Dokkum & Franx (1996) first noted, it is difficult to reconcile the high 
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Figure 3. The FP redshift estimates compared to the spectroscopic 
redshifts. The left panel shows the results for the lens galaxies, and 
the right for the cluster galaxies. For the lens galaxies, solid (open) 
points are used for lenses where the source redshift is known (unknown) . 
The dashed lines are offset by \Az\ = 0.1 to illustrate the desired 
accuracy. We have only H-band data for the labeled lenses with large 
redshift uncertainties (RXJ 0911+0551, MG 2016+112, B 2045+265 
and B 2319+052). 

Oo models with a reasonable star formation epoch (zf 5s 10). For the low S7o 
models, a coherent FP matching the local observations is possible for formation 
epochs Zf J> 2. The scatter about the FP is comparable to that for the local FP, 
so the lens galaxies cannot have an enormous range of formation epochs. 

Since we are now confident that the lenses lie on the FP, Figure 2 shows 
the evolution of the E+K corrections with redshift for both the lens and cluster 
samples in the f^o = 0-3 flat cosmological model. The results for an Qq = 0.3 
open model are similar. The dominant term depends on the wavelength. In the 
V band, the K-correction dominates and the galaxies rapidly fade, in the I band, 
the E and K corrections balance, and in the H-band, the evolution dominates 
and the galaxies steadily become brighter. At all bands, the galaxies become 
steadily brighter than predicted by a no evolution model. 

It is difficult to distinguish the evolution of the lens and cluster galaxies, 
and both samples favor star formation epochs of Zf ~ 2-3. This strongly con- 
tradicts the predictions of Kauffmann & Chariot (1998), where most of the field 
early- type galaxies form their stars at zj <J 1. Sommerville & Primack (1998) 
have argued more generally that the earlier semi-analytic models systematically 
underestimated the epoch of star formation through their choice of star forma- 
tion mechanisms. In their defense, Kauffmann & Chariot (1998) would argue 
that morphologically selected early-type galaxies look old at all redshifts be- 
cause the mergers required to produce the early-type morphologies do not occur 
when the stellar populations are young. This requires rapid number evolution 
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in the early- type galaxy population to z = 1, which is probably inconsistent 
with observations (Lilly et al. 1995, Schade et al. 1999). Moreover, the lens 
galaxies were selected based on mass, not morphology, and while we excluded 
a few late-type galaxies, there is no significant population of blue, star-forming 
lenses which we have dropped from the analysis. 

4. Photometric Redshifts 

Redshift incompleteness is the bane of many astrophysical applications of grav- 
itational lenses, even though there is no technical barrier to measuring most of 
the missing redshifts in the age of the 8 meter telescope. This is unfortunate, 
because you obtain more cosmological information with the measurement of a 
lens redshift than any other single object. For example, the distribution of lens 
redshifts is a powerful cosmological test (Kochanek 1992), but it requires high 
redshift completeness because it is very difficult to make a statistical model for 
incomplete redshift samples. The difficulty in measuring lens redshifts is further 
evidence that the stellar populations of lens galaxies are old. If the typical lens 
had any ongoing star formation, the redshift measurements would be far easier 
due emission lines usually associated with star formation. 

Photometric redshifts should work well for lens galaxies because of the re- 
stricted galaxy type distribution (mostly early-type galaxies, and no funny dwarf 
galaxies) and because we also know the galaxy masses. We have developed a 
redshift estimation method based on the FP. Placing a galaxy on the FP can 
accurately estimate a galaxy redshift with no color information if the lens has 
z\ £5 0.5 or we have optical photometry. If we have color information, partic- 
ularly a color bracketing the 4000A spectral break, then the color information 
constrains the redshift more strongly than the FP. Figure 3 tests the accuracy 
of the method by comparing the photometric redshift estimates to the known 
spectroscopic redshifts for the lens and cluster galaxies. Keep in mind that these 
estimates are made with 1-3 filter photometry. 

For the lens galaxies, all the poor estimates with large redshift uncertainties 
are lenses at z\ > 0.5 for which we have only infrared photometry. With the 
addition of any optical data the uncertainties collapse and the prediction usually 
matches the measured redshift. In a few cases (Q0142-100 and B1608+656 in 
Fig. 3), we have accurate 3-filter photometry and small redshift uncertainties, 
but the predicted redshift is incorrect. For B1608+656, the problem is that the 
galaxy is an E+A galaxy with more recent star formation (Myers et al. 1995). 
The average accuracies, {zfp — zt rU e), are —0.03 ± 0.10 and —0.02 ± 0.07 for the 
lens and cluster galaxies respectively. Much of the scatter for the lens galaxies 
is due to the systems with only H-band data. 

5. The Mass and Luminosity Function of Lens Galaxies 

The lens sample is large enough to begin studying the luminosity and mass 
functions of the lens galaxies. The average image separation of the lenses is 
already the most accurate way of estimating the mean velocity dispersion of 
massive galaxies, which makes matching the image separations a critical compo- 
nent of any estimate of the cosmological model using gravitational lens statistics. 
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Figure 4. The luminosity (left) and velocity (right) functions of lens 
galaxies compared to local clusters. The heavy solid (light dashed) his- 
tograms show the observed distributions of the lens (cluster) galaxies, 
and the heavy solid (light dashed) curves show the predicted distribu- 
tions for the lens (cluster) galaxies. The sharp feature in the predicted 
luminosity function for the cluster galaxies is created by our ad hoc 
model for the magnitude limits in the JFK sample of cluster galaxies. 



Because the lenses are a mass-weighted galaxy sample, they should be more lu- 
minous and more massive than a magnitude limited sample of cluster galaxies. 
We can make a preliminary comparison by assuming the standard model we 
have used for gravitational lens statistics (Kochanek 1996, Falco et al. 1998). 

We start with a standard Schechter luminosity function to count the galax- 
ies, 

dn/dL = (n*/L*)(L/L*) a eM-L/L*), (6) 

combined with a Faber-Jackson law to relate luminosity and velocity dispersion, 



L/L* = (a/a*)"' 



(7) 



where the velocity dispersion is related to the image separation through eqn. (5). 
The multiple imaging cross section, oc A# 2 oc a 4 oc L 4 ' 7 , is used to compute the 
probability that a galaxy acts as a lens. We fix the parameters to the standard 
model we have used to study lens statistics (7 = 4, a = —1, cr* = 225 km s _1 , 
and B* = —19.9 + 5 log h mag). The galaxy density n* will be nearly degenerate 
with the cosmological model, so we will not examine it here. 

Both the lens and cluster galaxy samples are affected by selection effects. 
For the local JFK cluster galaxy sample, it is the magnitude limit of the sample, 
which we have treated only approximately. For the lens galaxies, the only im- 
portant selection effect is the finite angular resolution of the surveys, which we 
include as a minimum detectable separation A9 m i n . The minimum separation 
leads to a minimum lens luminosity L m i n /L* = (A6 m i n /A9*)' y > 2 and a selection 



function S(x = L/L r , 



(10 — 15x + 6x 2 ). Combining all these effects and 



using our standard parameter values, the distribution of all early-type galaxies 
is 

dn/dM oc dn/d log a oc exp(— L/L*) L > L rn i n (8) 

where L m i n is set by the survey magnitude limit, and the distribution of lens 
galaxies is 

dn/dM oc dn/dloga oc (L/L*)exp(—L/L*)S(L/L min ) (9) 

where L m i n is set by the survey resolution. For the local cluster sample, the 
distribution of galaxies peaks at the luminosity L m i n corresponding to the sample 
magnitude limit. The lens galaxy sample peaks at L ~ L*, while the selection 
function matters only near L m i n ~ 0.05L*. In fact, our model for the selection 
function of the lens galaxies is better than that for the JFK cluster galaxies. 
The results are shown in Figure 4, where we present histograms of the observed 
lens and JFK cluster galaxy distributions as compared to our simple model. 
As expected, the lens galaxies are both more luminous and more massive than 
typical early- type galaxies. More importantly, the differences closely match the 
predictions of our standard model. 

6. Future Growth 

This is only an interim report on studies of the distribution and evolution of 
galaxies using gravitational lenses, since we have used only half the known lenses 
and the number of lenses is still growing rapidly. Our progress is limited by two 
major problems. First, despite our best efforts, only about 70% of the lenses 
have archival or scheduled HST images good enough to do surface photometry 
of the lens galaxy even under a generous definition of "good enough." Ac- 
curate measurements of galaxy evolution, mass-to-light ratios, luminosity and 
mass functions require higher precision data than the fractional-orbit snapshots 
commonly used to confirm lens candidates. Host galaxy contamination, while 
scientifically valuable (see Bernstein et al. and Rix et al. in these proceedings), 
is also a serious problem for many radio lenses. Second, the high level of redshift 
incompleteness is a major barrier to using gravitational lenses as astrophysical 
tools. Even so, the sample of 30 early-type galaxies we used to study galaxy 
evolution is 5 times larger than the biggest published sample of field early-type 
galaxies with kinematic measurements at comparable redshifts. We hope to 
solve the first problem with the continuation of the CASTLES project into HST 
Cycle 9, and the missing redshift problem should soon be solved by the advent 
of so many new large telescopes. 
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